Back

Journal of Bioinformatics and Systems Biology

Fortune Journals

Preprints posted in the last 7 days, ranked by how well they match Journal of Bioinformatics and Systems Biology's content profile, based on 15 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit.

1
Abundant Glomerular Neutrophil Extracellular Traps in C3 Glomerulopathy

O'Sullivan, K.; khandelwal, p.; Walker, P. D.; hickey, m.; Licht, C.

2026-09-01 immunology 10.64898/2026.08.27.747386 medRxiv
Top 0.3%
0.6%
Show abstract

Introduction: C3 glomerulopathy (C3G) is driven by fluid-phase alternative complement pathway dysregulation, with emerging evidence linking glomerular neutrophil infiltration to disease severity. Neutrophil extracellular traps (NETs) are implicated in other forms of glomerulonephritis. However, their participation in the pathogenesis of C3G remains undefined. Methods: Kidney biopsies from 33 patients with C3G (15 with dense deposit disease [DDD] and 18 with C3 glomerulonephritis [C3GN]) were compared with 15 anti-neutrophil cytoplasmic antibody associated vasculitis (AAV) biopsies as a neutrophil-rich disease control in this retrospective cross-sectional study. Glomerular neutrophils and NETs were identified using immunofluorescence, staining for myeloperoxidase, citrullinated histone H3, peptidyl arginine deiminase-4, and DNA. Supervised machine learning was used to quantify glomerular NET formation, and the data were correlated with kidney function at time of biopsy using linear regression. Results: Intraglomerular NETs were abundant and detected in the majority of glomeruli in C3G biopsies. Compared with AAV, C3G showed a significantly higher fraction of neutrophils forming NETs, despite similar neutrophil counts per glomerulus. NET abundance was similar in DDD and C3GN. In exploratory analyses, a greater proportion of glomeruli containing NETs was associated with lower kidney function (estimated glomerular filtration rate) at biopsy, and this association remained significant after adjustment for age, C3G subtype, and interstitial fibrosis. Conclusions: These observations demonstrate that intraglomerular NETs are a common and prominent observation in C3G and are associated with reduced kidney function at biopsy. These findings raise the possibility that NET deposition in glomeruli is a previously unrecognized driver of glomerular injury in C3G.

2
CyChat: a conversational Cytoscape app for no-code, reproducible network analysis

Liebold, J.; Stahl, M.; Schulze, J.-O.; Razavi, M. M.; Bader, G. B.; Kurtz, S.; Baumbach, J.

2026-09-01 bioinformatics 10.64898/2026.08.28.747833 medRxiv
Top 0.4%
0.6%
Show abstract

Network-based analyses of molecular interactions are useful for interpreting high-throughput omics data and identifying therapeutic targets. Cytoscape is the standard platform for these tasks, but users face a trade-off between accessible graphical workflows that are difficult to document and reproducible automation in Python or R that requires programming expertise. General-purpose coding assistants can generate Cytoscape Automation scripts, but remain external to Cytoscape. We present CyChat, a Cytoscape Desktop app that integrates a chat interface and a large language model (LLM) agent into the application. CyChat translates natural language into executable Cytoscape Automation workflows, runs generated Python code, and exports chat sessions with executed code as standalone Jupyter notebooks. To reduce setup barriers, CyChat includes an embedded Python runtime and supports both cloud-based and locally hosted LLMs. CyChat was evaluated across ten Cytoscape workflows using seven LLM providers, each represented by one LLM. The strongest configuration achieves a pass rate above 99%. In a qualitative evaluation based on a published network visualization, CyChat completes the task in 1.5-5 minutes, compared with 15-20 minutes for manual GUI workflows by computational biologists. CyChat is available through the Cytoscape App Store at https://apps.cytoscape.org/apps/cychat.

3
scPyviewer: a Python-native interactive viewer from AnnData single-cell data

Xuan, H.; Huang, Y.; Bian, J.; Liu, X.

2026-08-31 bioinformatics 10.64898/2026.08.26.747418 medRxiv
Top 1%
0.3%
Show abstract

Motivation: Interactive tools that let non-programmers explore an analyzed single-cell dataset, its embeddings, gene expression, cell metadata, and marker genes, have become standard laboratory infrastructure. Every actively maintained tool in this space (ShinyCell, ScRDAVis, sCIRCLE, scViewer) is built on R Shiny and requires a Seurat object as input. Laboratories whose primary analysis pipeline is Python/scanpy, the dominant framework for single-cell RNA-seq, spatial, and multi-omic analysis, therefore have no lightweight, language-native option that pairs a shareable web-based viewer with a scriptable Python API: sharing a scanpy result means either exporting to Seurat first or handing over a notebook that only a programmer can run. Results: We present scPyviewer, a web-based viewer that ingests AnnData objects directly and reproduces the core interaction patterns of the incumbent R Shiny tools without leaving the Python stack. In a feature-parity audit against three actively maintained R Shiny incumbents, scPyviewer matches or exceeds every baseline capability (7/7); among these, it uniquely offers native AnnData ingestion with no Seurat conversion, and cross-dataset comparison over shared genes and matched cell-type composition. Benchmarked head-to-head against the R/Seurat rendering substrate the incumbents are built on, identical operations, identical data, across three datasets spanning 22,315 to roughly 313,000 cells, scPyviewer renders every core view faster at every scale tested (up to 3.6x on a single view) and at a fraction of the memory (5.2x lower on the smallest dataset). At the largest scale tested, the gap becomes categorical rather than incremental: scPyviewer completes every view on a 313,000-cell dataset while the Seurat substrate exhausts an 8 GB memory budget and fails outright. Beyond the interactive app, scPyviewer installs via pip or conda and exposes a public Python API that returns Matplotlib figures and pandas tables for scripted, publication-ready output. Availability and implementation: scPyviewer is implemented in Python 3.11 (scanpy 1.11.5, anndata 0.12.19, streamlit 1.59.2, plotly 6.9.0) and distributed with a one-command reproduction interface that installs pinned dependencies, regenerates the benchmark and all figures, and launches the interactive app. Source code is available at https://github.com/xuan13hao/scPyviewer.git.

4
Automatic bioinformatic software named entity recognition from literature

Xuan, H.; Pasupuleti, R.; Liu, B.; Sun, H.; Zhang, J.; Yao, Z.; Zhong, C.

2026-09-01 bioinformatics 10.64898/2026.08.26.731133 medRxiv
Top 1%
0.3%
Show abstract

Bioinformatics software and databases are essential components of modern life science research, yet their mentions in the scientific literature are often inconsistent and difficult to systematically identify at scale. The lack of a comprehensive and up-to-date catalog of bioinformatics resources hinders efforts toward automated biomedical knowledge extraction and streamlined data analysis. Here we present SNAIL, a hybrid named entity recognition framework designed to automatically identify bioinformatics software and database (SW/DB) names from biomedical texts. SNAIL integrates complementary lexical and semantic modeling strategies. The lexical component captures orthographic patterns and contextual cues characteristic of SW/DB names, while the semantic component leverages contextual embeddings generated by transformer-based language models such as SciBERT, combined with an explicit token-masking strategy to enhance entity-focused representations. A large training corpus was constructed automatically through a hybrid pipeline that integrates citation-hinted extraction with large language model-assisted distillation. Evaluation on two independent benchmark datasets and real-world research articles demonstrates that SNAIL substantially outperforms existing approaches, including domain-specific methods such as bioNerDS2 and general-purpose large language models such as ChatGPT, Gemini, Grok and Claude. Applying SNAIL to large-scale literature analysis further reveals distinct journal-level preferences across bioinformatics subfields. These results demonstrate that SNAIL provides an accurate and scalable solution for identifying bioinformatics resources in scientific texts and enables systematic meta-analysis of tool usage and research trends.

5
A Randomized Non-Inferiority Trial of an eHealth Delivery Alternative for Cancer Genetic Testing for Hereditary Cancer (eREACH2)

Lee, K. T.; Egleston, B.; Fetzer, D.; Domchek, S. M.; Fleisher, L.; Wen, K.-Y.; Wagner, L.; Roberts, S.; Howe, S.; Cacioppo, C.; Christiansen, J.; Karpink, K.; Selmani, E.; Mastaglio, E.; Weinberg, M.; Wood, E. M.; Feng, J.; John, S.; Schweickert, K.; Mcleod, B.; Bradbury, A. R.

2026-09-03 genetic and genomic medicine 10.64898/2026.09.01.26361920 medRxiv
Top 2%
0.2%
Show abstract

Background: Many at-risk patients lack access to genetic services due to a genetic counselor (GC) workforce shortage. Little is known about how digital alternatives impact patients with and without cancer who meet criteria for genetic testing. Methods: eREACH2 is a randomized 4-arm non-inferiority trial where pre-test (visit 1) and/or return of results (visit 2) GC counseling was replaced with a patient-centered digital intervention. Arms include: A (GC/GC), B (GC/digital), C (digital/GC) and D (digital/digital). Primary outcomes were non-inferiority in uptake of genetic services and change in genetic knowledge and general anxiety from baseline to post-disclosure of results (T0-T2). Secondary cognitive and affective outcomes were assessed using non-inferiority ANOVAs and equivalency chi-squared tests in intention-to-treat and per-protocol analyses. Findings: 773 participants were recruited nationwide; 46.6% from rural areas. Mean age was 51 years (range 20-87), 13% male, 12% non-white, 29% had less than a college education, and 33% had a personal history of cancer. 584 (76%) patients completed testing (14% had a positive result, 16% had a VUS). In the primary ITT analyses, we met the non-inferiority for uptake of genetic services and anxiety, but results were inconclusive for knowledge. Secondary outcomes were heterogeneous across arms. Arm C demonstrated consistently favorable effects, while Arms B and D showed less favorable outcomes in select domains (e.g. satisfaction and MICRA). Patients who received positive or VUS results via digital disclosure had significantly higher MICRA scores - indicating greater negative response to testing. Interpretation: In this large, randomized trial of patients with and without cancer, the eREACH intervention was effective for pre-test counseling, but inconclusive for digital disclosure of results. Exploratory analyses suggest that digital delivery could be a reasonable alternative for individuals receiving negative results, while those receiving positive or VUS results may derive some short-term psychosocial benefit from GC disclosure.

6
Rclade: automated taxonomic collapsing and geological-timescale annotation of time-calibrated phylogenetic trees in R

Zeng, Z.; Wang, Y.

2026-09-01 bioinformatics 10.64898/2026.08.27.747462 medRxiv
Top 2%
0.2%
Show abstract

Background: Reproducible taxonomic collapsing and geological-timescale annotation of time-calibrated phylogenetic trees in R often require coordination among several packages and repeated code for label parsing, clade validation, plotting, and export. Workflow-managed analyses additionally benefit from non-interactive configuration, predictable diagnostics, and machine-readable exit status. Results: We present Rclade, an R package that consolidates the multi-package coordination required for taxonomic collapsing into a streamlined, single-function interface. Rclade provides (1) custom ggproto objects (GeomPolygonStraight/GeomSegmentStraight) that bypass coord_munch() interpolation to achieve straight-edge rendering of collapsed triangles in circular layouts; (2) automatic detection and parsing of four taxonomic-label formats (GTDB, Silva, NCBI, embedded) plus user-supplied custom regex, with explicit input-validation contracts and parsing-accuracy evaluation on real and derived test sets; and (3) workflow embeddability through YAML configuration, library-mode APIs, and standard Unix exit codes. Benchmarks on synthetic and real datasets (200-10,000 synthetic tips and real reference trees up to 10,122 tips; 5 replicates at every scale under a unified fully rendered measurement protocol) show that the full-pipeline overhead is modest for interactive use (median {approx}0.87 s in-session rendering and {approx}8.4 s process-level wall-clock at 10,000 tips). Conclusions: Rclade is a convenience layer over the ggtree/deeptime ecosystem that reduces boilerplate while adding targeted technical improvements for circular-layout rendering and format heterogeneity management.

7
From Prompt to Provenance: BloClaw, a Capability-Gated AI4S Workstation for Auditable Computational Biology

qin, y.; Pang, J.; Zhang, X.

2026-09-01 bioinformatics 10.64898/2026.08.26.747436 medRxiv
Top 2%
0.2%
Show abstract

Scientific agents can produce plausible answers while remaining unable to establish whether the computation behind an answer is executable, recoverable, or reproducible. We present BloClaw, an AI4S workstation built around a simple principle: a scientific agent should know what it can do, show how it did it, and state what remains unvalidated. Each capability declares an execution state, input constraints, dependencies, expected outputs, and scientific limitations. Natural-language requests are translated into structured tasks, validated against this registry, executed through scientific tools, and recorded in a provenance-aware Living Lab Notebook. The system is designed to detect invalid inputs, failed tool calls, missing dependencies, and remote timeouts, and to route them to repair, retry, or escalation. The implemented and tested scope comprises RDKit-based molecular property and rule screening, protein structure analysis, docking-pose inspection, 3D visualization, and structured reporting. We demonstrate the workflow on a PubChem-retrieved osimertinib structure and a supplied 6LU7 docking artifact: the former yields deterministic descriptors (molecular weight 499.619 Da, cLogP 4.5098, TPSA 87.55 A^2), while the latter contains 2,387 protein ATOM records, 309 residues, and nine pose records. These examples are workflow demonstrations, not efficacy or affinity studies. Beyond retrospective prediction, the manuscript specifies a prior-minimized constructive mode in which a desired function is compiled into explicit physical, chemical, and systems constraints, candidate mechanisms are simulated, and observations are reintroduced for calibration and falsification; this is a proposed extension rather than a result of the present case studies. We describe an evaluation protocol that compares BloClaw with a standard single-agent workflow and fixed-script execution using task completion, scientific correctness, recovery success, provenance completeness, reproducibility, human review time, latency, and cost. This manuscript reports the system design, verified capability boundary, deterministic software artifacts, and a reproducible evaluation protocol; it does not claim benchmark improvements before those experiments are run. BloClaw is an execution and accountability layer for AI-assisted research, complementing expert review and experimental validation rather than replacing them.

8
MetaDome 2027: a comprehensively updated resource for aggregating missense variant evidence across homologous human protein domains

Wiel, L.; Ferraro, F.; Yu, J.; Zhen, J.; Nachun, D.; Mendez, R.; Reuter, C. M.; Cui, J. L.; Bonner, D. E.; Carter, J. N.; Marwaha, S.; van de Vorst, M.; Emami, S.; Kravets, E.; Neu, M. B.; van Ham, T. W.; Kleefstra, T.; Ashley, E. A.; Bernstein, J. A.; Montgomery, S. B.; Gilissen, C.; Wheeler, M. T.

2026-08-31 bioinformatics 10.64898/2026.08.26.747388 medRxiv
Top 2%
0.2%
Show abstract

The interpretation of missense variants remains a major challenge in clinical genetics. "Meta-domains" aggregate population and pathogenic variation across homologous Pfam domain instances in the human proteome, providing per-residue context for interpreting variants of uncertain significance (VUS). Our 2019 implementation, MetaDome, is widely used and named in clinical variant-classification guidelines. Here we present the MetaDome 2027 update, featuring a comprehensively updated dataset and GRCh38 support. The redesigned pipeline enables incremental updates of GENCODE, UniProtKB/Swiss-Prot, Pfam, gnomAD, and ClinVar while maintaining 100% sequence-identity gene-to-protein mapping. Annotated Pfam domain instances grew 14.9% from 71,419 to 82,069 and meta-domain-eligible Pfam families ([≥]2 human occurrences) by 73.3% from 3,334 to 5,778; Pfam domains are annotated to 92% of human proteins. Approximately 43% of mapped protein-coding nucleotides (14.3 million in GRCh38, 13.8 million in GRCh37) are in a meta-domain; in GRCh38 67.9% (37,692 of 55,548) of pathogenic or likely pathogenic ClinVar missense variants fall at such a position. We show how MetaDome helped reclassify a de novo missense VUS in RALA and identify 52,463 ClinVar missense VUS for which meta-domains supply otherwise unavailable pathogenic evidence. MetaDome is freely available at www.metadome.app.

9
Genetic Architecture and Sample Size Impact Relative Performance of Nonlinear Machine Learning and Standard Polygenic Risk Scores

Zhu, J.; Baousi, A.; Morris, A. P.; Guo, H.

2026-09-03 genetic and genomic medicine 10.64898/2026.08.29.26361109 medRxiv
Top 2%
0.2%
Show abstract

Standard polygenic risk scores (PRSs) are constructed based on additive genome-wide association study (GWAS) summary statistics. Nonlinear machine learning methods have been increasingly applied to construct PRSs directly from individual-level data, with the aim of improving predictive performance over standard PRSs through their ability to model non-additive genetic effects. However, their superiority across studies has been inconsistent, and the conditions under which they provide meaningful improvements remain unclear. We combined theoretical analysis, simulations and a real-world application to investigate when two widely used nonlinear machine learning methods, random forest and XGBoost, outperform standard PRSs. Theoretical analysis showed that standard PRSs can implicitly capture part of the genetic variance attributable to nonadditive genetic effects through their contributions to marginal SNP effects, thereby losing less information than commonly assumed. Although nonlinear models have a higher theoretical potential, their greater flexibility incurs a bias-variance trade-off that can limit predictive gains at finite sample sizes. Simulations showed that XGBoost outperformed the standard PRS only when the genetic architecture involves a sufficiently large proportion of interaction genetic variance concentrated across relatively few interaction effects and large training samples were available. Random forest consistently underperformed the standard PRS. In an application to ischemic heart disease prediction using UK Biobank data, XGBoost showed no meaningful improvement in predictive performance over the standard PRS, whereas random forest again performed worse. Together, these findings suggest that nonlinear machine learning do not uniformly outperform standard PRSs; rather, their relative performance depends jointly on genetic architecture and training sample size. Our study helps to reconcile the inconsistent results reported across previous studies and provides a framework for identifying settings in which more complex PRS models are likely to be beneficial.

10
PhageTransformer - scalable and accurate host assignments for bacteriophages

Siemers, M.; Lopez, J. L.; Dutilh, B. E.

2026-08-30 bioinformatics 10.64898/2026.08.29.748026 medRxiv
Top 2%
0.2%
Show abstract

Bacteriophages can only be understood through their interactions with bacterial hosts. As environmental sequencing efforts expanded, the number of available phage genome sequences has exploded, yet the vast majority of these sequences lack host information. Predicting the host of a newly observed phage is therefore a key challenge in virology. Several computational tools can predict phage-host relationships from genomic data, but they share notable limitations: (1) the number of different hosts that can be predicted remains relatively restricted; (2) tools tend to assign confident host predictions to non-viral input sequences; and (3) most tools have a trade-off between accuracy and speed. Here we present PhageTransformer (PT), a deep learning model for phage-host prediction that addresses these limitations. We benchmark PT against existing tools on 3,881 independent phage-host pairs from GenBank and public HiC data, and demonstrate that it achieves competitive or superior prediction accuracy at greatly reduced runtime.

11
Relation of Self-Reported Race and Genetic Ancestry to Hypertension Prevalence Among Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos

Montanez-Valverde, R. A.; Kim, V.; Duran-Luciano, P.; Yuan, Y.; Sofer, T.; Kaplan, R. C.; Gallo, L. C.; Talavera, G. A.; Perreira, K. M.; Daviglus, M. L.; Rosas, S. E.; Llabre, M. M.; Elfassy, T.; Li, X.; Isasi, C. R.; Rodriguez, C. J.

2026-09-03 genetic and genomic medicine 10.64898/2026.09.01.26361995 medRxiv
Top 2%
0.2%
Show abstract

Background. The imprecision of current metrics to capture the complex genetic admixture and racial identity among Hispanic/Latino individuals in the United States [US] is a concern. We examined the relationship of self-reported race and genetic ancestry with hypertension [HTN] among Hispanics/Latinos. Methods. Cross-sectional study of the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), including 10,586 Hispanic/Latino unrelated adults. Genetic ancestry: West African [AA], Amerindian [AI], and European [EA]. Self-reported race: White, Black, Native American, or Multiple/Missing (More than one race or Unknown/Not reported/Refused). HTN: systolic (SBP) [≥]130 mmHg, diastolic blood pressure (DBP) [≥]80 mmHg, and/or use of HTN medications. Age- and sex adjusted models were used. Results. Self-reported race was White (38{middle dot}6%), Black (3{middle dot}6%), Native American (4{middle dot}1%), and Multiple/Missing (53{middle dot}7%), with Unknown/Not reported/Refused representing 32{middle dot}7%. Black and White Hispanics/Latinos had the greatest AA (55{middle dot}7%) and EA (69{middle dot}3%) ancestries, respectively. Each 10% AA increase was associated with OR 1{middle dot}15, SBP beta +0{middle dot}9 mmHg, and DBP beta +0{middle dot}7 mmHg. Conversely, each 10% AI increase was associated with OR 0{middle dot}83, SBP beta -0{middle dot}4 mmHg, and DBP beta -0{middle dot}6 mmHg. HTN prevalence was highest among those with Black race or in the highest AA quantile (45{middle dot}6% and 48{middle dot}0%, respectively), and lowest among those with Native American race or in the highest AI quantile (37{middle dot}6% and 26{middle dot}7%, respectively). Conclusion. One-third of Hispanics/Latinos did not self-report race. Black or White self-reporting race did somewhat relate to AA or EA ancestry, respectively. HTN profiles were related to self-reported race and genetic ancestry in this admixed population.

12
Mural-VISTA: A tool for mural cell-vessel interaction assessment and multiscale single-cell topo-morphological analysis

Zeng, H.; Hu, M.; Phng, L.-K.; Matsunaga, Y. T.

2026-09-01 bioinformatics 10.64898/2026.08.27.747487 medRxiv
Top 2%
0.1%
Show abstract

Three-dimensional (3D) mural cell morphology is heterogeneous and coupled to vessel geometry, however, measurements from two-dimensional (2D) maximum intensity projections (MIP) obscure overlapping processes and cell-vessel contacts. Accordingly, we developed Mural-VISTA, a semi-automated Python workflow for mural cell-vessel interaction and single-cell topo-morphology analysis of reconstructed surface meshes. This workflow integrates mesh pretreatment, interactive centerline extraction, hierarchical segmentation of cell soma, main axis and secondary processes (branches), and extraction of 36 multiscale (cell process segment level, process level, and whole cell level) topo-morphological and vessel-referenced metrics. Mural-VISTA identified morphological changes in pericytes and vascular smooth muscle cells (vSMCs) with altered RhoA activity. Constitutive active RhoA (RhoA CA) over-expression reduced branch complexity and increased process alignment in both cell types, while increased whole-cell and branch solidity only in vSMCs. Dominant negative RhoA (RhoA DN) over-expression increased branch abundance and reduced branch solidity in pericytes but not vSMCs, suggesting cell-type specific effect of reduced RhoA activity. In conclusion, Mural-VISTA enables quantitative 3D profiling of mural cell architecture and its spatial relationship with the vessel.

13
Network-based meta-analysis maps stage-dependent molecular programs in MASLD through MASLD-META NETWORK application

Kumak, E.; Darde, T.; Konu, O.

2026-08-31 bioinformatics 10.64898/2026.08.26.747338 medRxiv
Top 2%
0.1%
Show abstract

Metabolic dysfunction-associated steatotic liver disease (MASLD), the leading cause of chronic liver pathologies worldwide, represents a growing clinical burden. Its diagnosis remains reliant on liver biopsy that limits early detection and the ability to capture molecular changes across disease progression. A systematic understanding of stage-dependent gene expression changes is essential to identify biomarkers and effectively characterize disease mechanisms. Therefore recent studies provided databases for searching genes as well as prediction of multi-gene signatures for disease progression. However, there is still a need for interactive and comprehensive meta-analysis of datasets of MASLD patients with available histological metadata. Herein, we performed a meta-analysis of RNA-seq datasets using NAFLD Activity Score (NAS; n = 897) and fibrosis stage (n = 856) upon conducting pairwise comparisons across histological stages and identified differentially expressed genes associated with disease progression. Most importantly, we provide our findings via a dedicated web server, the MASLD-META NETWORK (https://masld.scilicium.com), enabling users to interactively explore meta-analysis results across diverse network modalities. In addition, we characterized gene expression dynamics across increasing disease stages to identify consistent progression-associated pathways using Louvain clustering. Network-based parameters such as centrality in combination with meta-analysis scores further highlighted central genes and pathways implicated in disease mechanisms. Accordingly, MASLD-META NETWORK enabled an integrative reassessment of recently published gene signatures, identifying COL1A1, COL3A1, THBS2, FBLN5, and PDGFA as the most central genes, and SULF2, MMP14, IL32, GPNMB, and COL3A1 as candidate markers of earlier transcriptional alterations. Network analysis of MASLD associated biological modules further identified LAMA2 and LAMA3 as previously unrecognized central candidate targets.

14
Dynamic Hierarchical Interleaved Bloom Filter: An Updatable Index for Large-Scale Fast Sequence Search

Seiler, E.; Willemsen, M.; Piro, V. C.; Reinert, K.

2026-08-30 bioinformatics 10.64898/2026.08.26.747224 medRxiv
Top 3%
0.1%
Show abstract

Motivation: A continued decrease in sequencing costs has facilitated the exponential increase in available sequencing data, with public databases like the European Nucleotide Archive (ENA) and Sequence Read Archive (SRA) reaching well in the order of petabases. This has been the incentive to develop more scalable tools for common bioinformatics tasks. One such task is the approximate searching of short sequence patterns like genes or reads in reference data sets. In recent years, a variety of indexing data structures have been proposed for searching large sequencing databases. The state-of-the-art index, the Hierarchical Interleaved Bloom Filter (HIBF) was first-in-class to index one million samples. To be useful for expanding repositories, it must be extended to support dynamic updates. Results: In this paper, we introduce a scalable and updatable sequence-search index by extending the HIBF with partial rebuilding to support efficient updates. We demonstrate the Dynamic HIBF's capacity for large-scale data by iteratively creating an index from over 100 TB of compressed reads across more than 39,000 full human RNA-Seq samples, updated in consecutive batches of 100. To benchmark against state-of-the-art tools, we evaluated incremental performance on a subset of 5,000 samples sub-sampled to 1% of their original read depth. In this comparative setting, the dynamic HIBF completed the sequential insertion of all 5,000 samples within 5 hours--24 to 65 times faster than competing methods and twice as fast as the static HIBF.

15
AmPair: automating housekeeping-gene primer design for species-level metataxonomics

Xu, X.; Yang, X.

2026-09-01 bioinformatics 10.64898/2026.08.25.746527 medRxiv
Top 3%
0.1%
Show abstract

Amplicon sequencing of the 16S rRNA gene is the most widely used approach for profiling bacterial communities, but its taxonomic resolution is typically limited to the genus level. Many species carry multiple divergent 16S rRNA alleles that overlap across species boundaries, an ambiguity that even full-length, long-read sequencing cannot fully resolve. Shotgun metagenomics achieves species-level resolution but remains costly, particularly when only a single genus is of interest. Amplicon sequencing of rapidly evolving, protein-coding housekeeping genes offers a cost-effective alternative, yet no tool exists to identify suitable primer sets for a given target taxon. Here we present AmPair, a Snakemake pipeline that, given a target genus and one or more candidate housekeeping genes, designs and ranks primer pairs binding conserved regions while flanking a variable region capable of species-level discrimination, and validates them in silico across all available genomes. Using the genus Bacillus and the housekeeping gene tuf as a case study, the primer set recommended by AmPair amplified 99% of 2,392 genomes; only 0.04% carried multiple alleles and none showed inter-species allele overlap, compared with 91.41% and 69.49%, respectively, for the standard 16S rRNA V1-V9 region. Applied to a Bacillus community profiled by Nanopore sequencing, the same primers resolved closely related species. AmPair thus offers a generalizable and accessible route to species-level community profiling.

16
A Curated Pharmacogenomic Allele Catalog for Sub-Saharan African Populations

SULAIMAN, M. A.; Oyeyemi, B. F.

2026-08-31 genetic and genomic medicine 10.64898/2026.08.25.26361354 medRxiv
Top 3%
0.1%
Show abstract

Sub-Saharan African populations carry pharmacogenomic alleles poorly represented in the European-derived reference panels underlying most clinical genotyping tools. We present a curated, machine-readable catalog of nine actionable alleles across six pharmacogenes (CYP2D6, CYP2B6, CYP2C9, CYP2C19, CYP3A5, NAT2) with African-specific frequency ranges, functional annotations, and evidence levels derived from reanalysis of 661 high-coverage whole-genome sequences across seven 1000 Genomes Project African populations. Direct comparison against PharmCAT v3.4.0 shows that CYP2D6 produces zero diplotype calls (0/661 samples callable) due to monomorphic reference positions absent from standard variant-only VCF output, a known limitation whose consequences for African allele carriers had not been reported. afripharmagen's reduced-position strategy identifies 243 CYP2D617 and 134 CYP2D629 carriers from the same input. For CYP2B6, CYP2C9, CYP2C19, and NAT2, both tools show concordance of 95-100%. Frequency gradients (CYP2B66: 30-50%; CYP2D617: 15-35% in West Africa; CYP3A5*1: 60-95%) translate directly into prescribing risk for efavirenz, tramadol, tacrolimus, and isoniazid. Pharmacogenomic decision support in African settings must incorporate population-specific allele definitions and input-format-aware strategies.

17
Development of iPSC-derived urothelial organoids towards investigating the effect of hormones on host-defense to urinary tract infections

Bindas, A.; Fang, Z.; Boekhorst, J.; Fernandes, A. M.; Wells, J.

2026-08-31 cell biology 10.64898/2026.08.29.747866 medRxiv
Top 4%
0.1%
Show abstract

Recurrent urinary tract infection represents a substantial unmet public health in women. Local administration of estradiol has been shown to reduce recurrence, however in vitro models of the female urinary tract remain limited and the mechanisms underlying the effects of estradiol are incompletely understood. Here, we describe a novel iPSC organoid differentiation protocol and its application to establish a multilayered transwell barrier culture model. Estradiol treatment resulted in reduced expression of innate antimicrobial peptides and cytokines, together with increased expression of demannosylation pathways. Treatment of transwell cultures with a combination of female sex hormones reduced endogenous CXCL8 signaling, independently of a 24-hour uropathogenic Escherichia coli (UPEC) challenge. To our knowledge, this is the first iPSC organoid-derived model of the urinary tract, which provides a platform for investigating interactions between the urothelium, urobiome and hormonal environment.

18
Design and Validation of New Primers for Specific and Sensitive Real-time PCR Detection and Quantification of Seven Botulinum Encoding Genes (Serotype A-G) of Clostridium botulinum

Phan, P.-L.; Chu, H.-A.; Le, T.-T.; Le, P.-A.; Nguyen, H.-L. T.; Tran, M.-N. T.; Nguyen, T.-T.; Pham, Y.; Phan, T.-N.

2026-09-01 molecular biology 10.64898/2026.08.21.746353 medRxiv
Top 4%
0.1%
Show abstract

Botulinum neurotoxins (BoNTs) comprise a highly diverse group of seven serotypes (from A-G) and over 40 subtypes worldwide. Previous primer- and probe-based nucleic acid amplification tests (NAATs) for detection of BoNT encoding genes are challenged by high levels of nucleotide polymorphism both across and within subtypes. In this study, multiple BoNT gene sequences were aligned to identify highly conserved regions for the design of new primers that enable the detection of all seven serotypes under the same conditions. Specific primer sets were designed and validated using in silico, conventional and real-time PCR with constructed plasmids carrying the target fragments and spiked food matrices. The established procedure achieved highly specific and sensitive detection of BoNT serotypes A-G with sensitivity of 10 copies/reaction and a total turnaround time of approximately 1.5 hours. The procedure also eliminated the carryover PCR product by using uracil-N-glycosylase in combination with dUTP in the assay reaction mix. This study provides an alternative NAAT with higher coverage and compliments the traditional mouse bioassays in enhancing global botulism surveillance capabilities.

19
Generation and characterization of a patient-specific human induced pluripotent stem cell line from a Skogholt syndrome patient (ASCFi003-A)

Przybyla, W.; Gupta, S.; Fjerdingstad, H. B.; Selnes, P.; Sharma, K.

2026-08-31 cell biology 10.64898/2026.08.29.747981 medRxiv
Top 4%
0.1%
Show abstract

We report the generation and characterization of a human induced pluripotent stem cell (iPSC) line derived from dermal fibroblasts of a patient with Skogholt disease, a rare maternally inherited neurodegenerative syndrome associated with choroid plexus dysfunction and impaired cerebrospinal fluid (CSF) homeostasis. Patient fibroblasts were reprogrammed using the non-integrating Repro-OSKGM kit. The resulting iPSC line exhibited typical pluripotent morphology, expressed canonical pluripotency markers, maintained a normal karyotype, retained the disease-associated genetic variant, was mycoplasma-free, and demonstrated trilineage differentiation potential. We also made choroid plexus (ChP) like organoids from the generated iPSCs. This patient-specific iPSC line provides a valuable resource for generating choroid plexus organoids and neurons to investigate disease mechanisms and develop therapeutic strategies.

20
Inhibition of JEV infection using β-Catenin specific inhibitor, iCRT-14

Datey, A.; Ghosh, S.; Chatterjee, S.; Bhowmick, B.; Ghatak, A.; Subudhi, B. B.; Chattopadhyay, S.

2026-08-31 molecular biology 10.64898/2026.08.29.747967 medRxiv
Top 4%
0.1%
Show abstract

The lack of effective anti-JEV therapy possesses significant challenge to control JEV. {beta}-catenin, a key mediator of Wnt signaling pathway regulates different viral replication and host immune responses. However, its role in JEV infection remains to be elucidated. Thus, the current study focused on evaluating iCRT-14, a specific {beta}-catenin inhibitor, against JEV. Treatment with iCRT-14 following JEV infection resulted efficient reduction in viral progeny release, viral RNA and protein levels in Huh7 and HEK293T cells. Further, active and total {beta}-catenin, Cyclin D-1 and GSK3-{beta}, the other key pathway players were also modulated in infected and inhibitor treated cells. Moreover, iCRT-14 showed an IC of 4.56 in Huh7 cell and maximal inhibition at the early stages of the JEV life cycle. Interestingly, the overexpression of {beta}-catenin in both the cells and siRNA-mediated {beta}-catenin knockdown (in Huh7 cells) significantly abrogated JEV replication, as evidenced by decreased viral titers, viral protein expression, and viral as well as total RNA levels. Moreover, the reduction in extracellular (84%) and intracellular (60%) viral titers following iCRT-14 treatment highlights its role in impairing JEV infection. Further, in silico molecular docking and co-immunoprecipitation studies demonstrated interactions between {beta}-catenin and the JEV NS5 and E proteins. Collectively, these findings suggest that optimum level of {beta}-catenin is required for efficient JEV infection, highlighting its potential as a target for designing host-directed control strategies to regulate viral infection.